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ABSTRACT 

Dwarf galaxies are widely believed to be among the best targets for indi- 
rect dark matter searches using high-energy gamma rays; and indeed gamma-ray 
emission from these objects has long been a subject of detailed study for ground- 
based atmospheric Cherenkov telescopes. Here, we update current exclusion 
limits obtained on the closest dwarf, the Sagittarius dwarf galaxy, in light of re- 
cent realistic dark matter halo models. The constraints on the velocity-weighted 
annihilation cross section of the dark matter particle are of a few 10 -23 cm 3 s~ 1 
in the TeV energy range for a 50 h exposure. The limits are extrapolated to the 
sensitivities of future Cherenkov Telescope Arrays. For 200 h of observation time, 
the sensitivity at 95% C.L. reaches 10 -25 cm 3 s _1 . Possible astrophysical back- 
grounds from gamma-ray sources dissembled in Sagittarius dwarf are studied. It 
is shown that with long-enough observation times, gamma-ray background from 
millisecond pulsars in a globular cluster contained within Sagittarius dwarf may 
limit the sensitivity to dark matter annihilations. 

Subject headings: Gamma-rays : observations - Dwarf Spheroidal galaxy, Dark 
Matter 



1. Introduction 

Dark matter (DM) plays a key role in the dynamics of a large class of astrophysical 
systems in the Universe. Though halos of dark matter are predicted to exist around all 
galaxies, dwarf spheroidal galaxies in particular are ideal targets for DM annihilation searches 
because: (i) their stellar dynamics show that they are among the most DM-dominated objects 
in the Universe; (ii) due to the lack of recent star formation activity, their environment is 
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relatively quiet in terms of background astrophysical gamma-ray emission ; (iii) many of 
them lie at distances below 100 kpc from the Galactic Center. 

The search for secondary gamma-rays from annihilations of dark matter particles is 
a powerful indirect detection technique because gamma-rays do not suffer from propaga- 
tion effects, the gamma-ray signal should be proportional to the square of the DM density, 
and characteris tic features such as lines or steps may be present in the energy spectrum at 
these energies (IBergstroml |2000| ; iBringmann et a/.ll2008i) . Imaging atmos pheri c Cherenkov 
telescopes (IACTs) such as HESS fao03j), MAGIC (Eofli) and VERITAS (12007h . are partic- 
ularly well suited to deep searches of targeted objects because of their large effective areas 
(~ 10 5 m 2 above 100 GeV). However, since IACTs are multipurpose astrophysical experi- 
ments and have a short duty cycle (~1000 hours/year), the observation time dedicated to 
these objects is typically limited to tens of hours per year. 

Since the flux of the expected gamma-ray signal is inversely proportional to the square of 
distance, one would expect the best dwarf spheroidal target to be the nearest one. However 
such dwarfs are also the closest to the Galactic Center and experience the tidal effect of 
the Milky Way. Recently, it has been shown t hat one could take advan tage of this effect to 
trace back the evolution history of the object (IPenarrubia et a/.ll2008a| ). During the orbital 
motion of a dwarf galaxy, multiple crossings of the dwarf galaxy through the galactic disc of 
the Milky Way give rise to the formation of tidal streams, a careful study of which allows to 
one infer the gravitational potential of the dwarf galaxy. 

In the case of the Sagittarius Dwa rf galaxy (SgrDw), the tidal streams have been de 
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) and have been used to derive the DM halo potential. Furthermore, 



measurements of stars within SgrDw and the luminosity of its core and surrounding debris 



allows the estimate of the D M content prior to tidal disruption (INiederste-Ostholt et al 



20101 ; IPenarrubia et aZ.I l2010h . Other peculiar features of SgrDw incl ude the presence of the 



M54 globular cluster coincident in position with its center of gra yitv (llbata et q/.lll994f) . and 



hints for the presence of a central Intermediate Mass Black Hole ( llbata et aZ.ll2Q09h (IMBH) . 
The latter point is supported by the observation of a deviation from a flat behavior in the 
surface brightness density profile towards the center of the object. 

Constraints on a DM annihilati on signal towards SgrDw, Canis Major, Sculptor an d 
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Carina have been reported by HESS (lAharonian et aZ.M2008l 
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Aleksic et aZ.I l201lh and towards Draco, Ursa Minor, Bootes 1 and Willman 1 by VERI- 
TAS ( lAcciari et a/.ll2010[ ). and towards Draco and Ursa Minor by Whipple flWood et al 
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20081 ). Because of its location in the Southern hemisphere, HESS is better suited for obser- 
vations of SgrDw with respect to other currently-operating IACTs. Observations with the 
Fermi Large Area Telescope (Fermi- L AT, a space-based telescope sensitive to gamma-rays 
between 20 MeV and 300 GeV), are also well-suited due to its large duty cycle and wide 
field-of-view, though the energy range probed is lower than that of IACTs. The Fermi-LAT 
collaboration put strong co nstraints in the GeV DM mass range on d warf spheroidal galaxy 
satellites of the Milky Way f Abdo et al. 2010a ; Ackermann et q/lboil ). However, their study 
is restricted to high galactic latitude (|b| > 30°) objects to avoid systematic contamination 
from galactic diffuse gamma-ray emission, and therefore no constraints on the measured flux 
in the direction of SgrDw have yet been published using Fermi-LAT data. 

In this paper, the current constraints on a DM annihilation signal toward s SgrDw are 



reass e ssed in light of more realistic DM halo models than previously used (jEvans et al 



2004 ; lAharonian et al\ 120081 ) . The s ensitivity of the future generation of IACTs, i.e. CTA 
(Cherenkov Telescope Array, |2010| ). is used to evaluate its potential for the detection of 
a DM annihilation signal. The CTA design-study sensitivity is used to investigate possible 
conventional gamma-ray emission, e.g. to the population of millisecond pulsars (MSP) in 
the globular cluster M54 at the center of SgrDw, or from the jet of a hypothetical central 
IMBH. It is shown that such standard astrophysical signals may limit the sensitivity to DM 
annihilations with CTA in case of long observation times, eventually requiring the modelling 
and subtraction of these astrophysical components. 

The paper is structured as follows: Section [2] is dedicated to the description of current 
and future instruments as well as the calculation of the sensitivity to DM signals. In Sec- 
tion [3j the modelling of the DM halo of SgrDw is described together with the astrophysical 
contribution to the DM flux. In the absence of an astrophysical gamma-ray background, 
exclusion limits on the velocity-weighted annihilation cross section of DM are derived in 
Section HJ Section deals with the estimate of the gamma-ray emission from the MSP 
population and the IMBH candidate of M54. Section 6 is devoted to the summary. 



2. Dark matter searches with IACTs 

2.1. Current and future instruments 

The present generation of IACTs (HESS, MAGIC and VERITAS) consists of multiple- 
telescope arrays detecting very high energy (VHE, E 7 > 100 GeV) gamma-rays. The stereo- 
scopic view of extensive air showers generated in the atmosphere by VHE gamma-rays al- 
lows these instruments to accurately reconstruct the direction and the energy of the primary 
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gamma-ray. The angular resolution reaches 0.1° per gamma-ray event and the point source 
sensitivity is about a f ew percent of the Crab Nebula flux above 100 GeV (see, for instance, 



Aharonian et al\ 120061 ). 



The plan for the next generation of IACTs, the Cherenkov Telescope Array (CTA, 120101 ). 
involves building two large arrays, one in each hemisphere, with an order of magnitude more 
telescopes than current instruments. This future instrument is expected to increase the flux 
sensitivity by a factor of 10 compared to current instruments, and enlarge the accessible 
energy range both towards the lower and higher energies. Based on the current CTA design 
study, a factor of about ten in effective area and a factor of two better in hadron rejection 
are expected. In this study, the estimated CTA ef fective area at the trigger level (before 
offline gamma-hadron separation) is extracted from |Paz Arribas (ESS). In order to mimic 
the effect of the analysis event selection, the effective area values for energies from ~100 GeV 
down to ~20 GeV are realistically lowered. The effective area then decreases from ~ 10 6 m 2 
at 200 GeV down to ~ 10 3 m 2 at about 20 GeV. 



2.2. Sensitivity calculation and background estimates 

The sensitivity for IACTs is calculated by comparing the number of events expected from 
an assumed gamma-ray emission scenario with the expected level of background events. In 
the case of DM searches, the assumed emission is from the annihilation of DM particles of 
mass m in the halo of the host galaxy, the differential gamma-ray flux of which is given by: 



d$(Aft,£ 7 ) 1 (av) dV 7 



x J(AQ)AQ, (1) 



d£ 7 8tt m? dE 1 

Particle Physics 

where (av) is the velocity- weighted annihilation cross-section and dN^/dE^ the photon spec- 
trum per annihilation. The astrophysical factor is defined as 

J(An) = -±- [ dtt [ p 2 [r(s)}ds. (2) 

When treating the self-annihilation of DM particles, this factor scales with the squared 
density of DM, p 2 , over the whole observation cone. The integral is then taken along the 
line of sight (LOS) and inside the solid angle Afl The solid angle is chosen as the angular 
resolution for point-like searches. The number of expected signal events can be calculated 
by: 

r°° d$ 

N 1 = T ohs A^E^—LdE^ (3) 
Jo dE 1 
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where T Q b s is the observation time, and A e g(E^) is the effective area of the detector as a 
function of the gamma-ray energy. In the case where the background is not measured exper- 
imentally, it can still be estimated assuming that the background consists of misidentified 
hadron showers. The estimate of the expected number of b ackground events in th e signal 
region can be determined using the following expression (see iBergstrom et q/.lll998l ): 
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where €h a d is the hadron detection efficiency To take into account the performance of 
the future IACTs the hadron rejection is taken at the level of 90%, which corresponds 
to ehad = 0. 1. This parametrisa tion gives remarkable agreement with CTA background 
simulations (IDi Pierro et q.Z.1 1201 lh . 



In case of no gamma-ray signal, a limit on the number of ga mma rays at 95% c onfi- 
dence level (C.L.), N* 5%aL -, can be calculated using the method of iRolke et al\ ( 120051 ). In 
what follows two cases are considered. In the case of current IACTs, the N^ 5%c - h - calcu- 
lation uses the numbers of gamma-ray and background events extracted from llh H.E.S.S. 
measurements ( lAharonian et aZ.ll2QQ8h . The projected _/V^ 5%C L ' for 50 h observation time is 
obtained by extrapolating both the numbers of gamma-ray and background events from llh 
to 50 h. In the case of 95% C.L. sensitivity calculations, J\f® 5%c - L - is calculated assuming the 
background-only hypothesis. For the H.E.S.S. sensitivity the number of background events 
is taken from the extrapolation at 50 h of observation. For the CTA sensitivity, the number 
of background events is calculated by integrating the background event flux given in Eq. (J3J) 
after multiplication by the effective area of the detector and the observation time. _/V^ 5%C L ' 
is then calculated using five off regions. 

Replacing Eq. ([I]) in Eq. 03]), the DM sensitivity can be then expressed in terms of the 
remaining particle physics parameters, (crv), m and dN/dE^. The 95% C.L. limit on the 
velocity-weighted annihilation cross section is given by the following expression: 
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3. Modelling the Sagittarius dwarf dark matter halo 



The Sagittarius dwarf (Sgrdw) is the only satellite 



evidence of o ngoing tidal mass stripping (llbata et al. 



tidal stream (IMateo et all Il998l ; iMajewski et al 



1999 



galax y in the MW that shows clear 



200 lh in the form of an associated 
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2004; Belokurov et al. 2006; Watkins et al 



2009) . This galax y is currently located at a close 
) . Indeed, it underwent its las t 



Mateo et al 



1998 



distance from the MW centre (~ 17 kpc; 
perigalacticon passage only 17 Myr ago f ILaw and Majewskill2010l ; IPenarrubia et a/.ll2009l ). 
which is a relatively short time compared with its internal dynamical time tdyn = R r J^ o ~ 47 



Myr, where R c is the galaxy core radius and ctq the central velocity dispersion (jMateolll998l ). 



The proximity of the Sgrdw to the Milky Way plus the fact that this galaxy is shed- 
ding stars to tides complicates its dynamical modelling in a number of ways. On the one 
hand, the distribution of dark matter and stars has been clearly altered from its original 
configuration by tidal mass strippin g. Given that the act ual amount of stars and dark 



matter in the tidal tails is unknown flNiederste-Ostholt et al.ll2010l ). the original mass, lu- 
minosity and size of the Sgrdw remain fairly uncertain quantities. On the other hand, the 
assumption of dynamical equilibrium may not be adequate, specially in the outskirts of the 
galaxy wh ere the population of u nbound stars may dominate in number over that of bound 
members ( iPenarrubia et a/.ll2009l ). 



These difficulties have not deterred a large body of theoretical work devoted to uncover 
the actual content and distribution of DM in the Sgrdw. To date these efforts have focused 
on (i) analytical models of the dynamical properties of the remnant core and (ii) N-body 
simulations that aim to reproduce the spatial and kinematical distribution of the tidal tails. 

The simplest analytical models assume dynamical equilibrium and adopt a cosmologically- 
motivated halo density profile to describe the kinematics of individual stars 

PnfwO) = , i Vl Ps , — , (6) 
{r/r s )(l + r/r s y 

where r s is a scale radius and p s is a characteristic density (Navarro, Frenk & White 1996, 
hereafter NFW). Note that this profile diverges at small radii as p oc r -1 , w hich is typically 
referred as a dark matter "cusp" . It was shown in IPenarrubia et all (|2008al ) that the tightly 
bound dark matter cusp is more resilient to disruption than th e mor e loosely bound stella r 
cored profile, which can be accurately described with a King (119661 ) profile ([Mated Il998l ). 
and that tidal stripping does not change the inner profile of DM haloes. 

Assuming that the external tidal field does not influence the kinematics of stars that 
locate the central regions of the dwarf, and ignoring the effects of tidal stripping on the 
outer (r ^> r s ) dark matter halo profile, one can use the Jeans equations to search the 
DM halo parameters that best fit the stellar central velocity dispersion for a observed King 
core radius of this object. The King-NFW degener acy gives rise to a family of NFW halo 
models which can reproduce the stellar dynamics (IPenarrubia et al\ l2008bl ). One way to 
break this degeneracy is using the relationship between th e virial mass and co ncentration 
found in cosmological N-body simulations (see for instance, iBullock et a/.ll200l[ ). Using this 
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procedure on the SDSS survey data provides a value of r s = 1.3 kpc. Considering the scatter 
on the relationship between virial mass and concentration, the la error on r s is found to 
be ~0.2 kpc. This correspond to the family of models with p s spanning from 7.5 x 10~ 3 to 
1.3 x 10~ 2 M Q pc~ 3 . In TabledJwe show the results of our fits together with the astrophysical 
factors J for different solid angles AQ. Taking into account the error on the halo profile 
parameters the value of the astrophysical f actor can vary by a factor of 2. Interestingly, an 
independent analysis by lLokas et al\ (120101 ) provides similar values for these parameters. In 
this case the astrophysical factors are found to be of a few higher than the ones presented 
here. 

However, numerical N-body models that aim to describe the observed structural and 
kinematical distributions of stars in the tidal tails as well as the remnant core provide a 
more consistent approach to the dynamical analysis of the Sgrdw. Yet, most of the existing 
N-body models of this galaxy assume for simplicity that dark matter and stars share the 
same spatial distribution (the so-called "mass-follows-light models"), an assu mption that is 



not su pported by detailed kinematic data of dwarf spheroidal galaxies (e.g. IWalker et al. 



(120091) ). The onl y excep tion to date corresponds to recent N-body models constructed by 
Penarrubia et all (120101 ). who explore the possibility that the Sgrdw may have originally 
been a rotating galaxy. In these models the galaxy is composed of an exponential stellar disk 
embedded in an extended DM halo. The DM density profile is taken as a cored isothermal 
(ISO) profile 

m h a exp[-(r/r cut ) 21 



Pisoir) = — T , (7) 

where is the halo mass, r c is the core radius and a ~ 1.156 (jPenarrubia et q/.ll2010l ). The 
DM halo mass can be estimated using the initial lumi nosity and a given mass-to-light ratio. 
Using the results from iNiederste-Ostholt et all (120101 ) the initial luminosity i s estimated t o 



be ~ 10 s L Q . Assuming a typical mass-to-light ratio for dwarf galaxies of 25 ( lMateolll998l ) 



the DM halo mass is found to be = 2.4 x 1O 9 M . To account for the initial tidal 
disruption of the SgrDw halo by the Milky Way, a truncation of the halo profile is imposed 
at r cut = 12 r c . The evolution of the SgrDw in the Milky Way potential is obtained via a 



N-bo dy model of SgrDw using the particle-mesh gravity code SUPERBOX (IFellhauer et al. 



20001 ) . The evolution code allows to recover the actual DM profile by using the constraint 
of the observed stellar distribution. The values of the parameters of the present ISO profile 
are given in Table [TJ 
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Exclusion limits on the dark matter annihilation cross section 



Theories beyond the Standard Model (SM) of particle physics propose several parti- 
cle DM candidates. For instance, some supersymmetric extensions of the SM predict a 
neutr alino as the lightest stable supersymm etric particle, which is a good candidate for 
DM (jJungman et all Il996l ; iBergstroml |2000| ). The para metrization of the neut ralino self- 
annihilation gamma-ray spectrum diV 7 /d.E 7 is taken from iBergstrom et al\ ( 119981 ) for a typ- 
ical neutralino annihilating into W and Z pairs. Fig. [T] shows the upper limits of current 
IACTs on (cry) as a funct ion of the DM mass m for AQ = 2 x 10~ 5 sr. Using the HESS 
upper limits published in lAharonian et all ( 120081 ) . the new upper limits are calculated for 
the NFW and ISO DM halo profiles of Section [3] and 11 h of observation time; the pro- 
jected upper limits for 50 h of observation time is also plotted. The limits are at the level of 
5 x 10 -23 cm 3 s _1 around 1 TeV for 50 h. The sensitivity of H.E.S.S. for 50 h observation time 
is also displayed. The sensitivity limits for CTA on (av) as a function of the DM mass m 
are presented in Fig. [3J for 50 h and 200 h observation times. The limits are calculated with 
AO = 2 x 10~ 6 sr for the NFW DM halo profile and AQ = 10" 3 sr for the ISO DM halo 
profile. The sensitivity limits at 95% C.L. reaches the level of 10 -25 cm 3 s _1 for DM masses 
of about 1 TeV in the case of the ISO DM halo profile. 

Two additional contributions to the overall gamma-ray flux that can modify the limits 
are considered: namely the Sommerfeld effect and Internal Bremsstrahlung (IB) from the 
DM annihilation. The Sommerfeld effect is a non-relativistic effect which arises when two 
DM particles interact in an attractive potential. When the relative velocity between the DM 
particles is s ufficiently low, the Som merfeld effect can substantially boost the annihilation 
cross section (ILattanzi and Silkll2009f ). since it is particularly effective in the very low- velocity 
regime. The actual velocity-weighted annihilation cross section of the neutralino can then 
be enhanced by a factor S defined as 

(av) = S(av) , (8) 

where the value of S depends on the mass and relative velocity of the DM particle. Assuming 
that the DM particles only annihilate to a W boson, the attractive potential created by the Z 
gauge boson through the weak force before annihilation would give rise to an enhancement. 

re stars, the 
). The value 



Assuming that the DM velocity dispersion inside the halo is the s ame as for t 



value of the DM velocity dispersion is fixed at 11 kms 1 f or SgrDw (IMateo 



1998 



of the enhancement is numerically calculated as done in ILattanzi and Silk! (120091 ) and then 
used to improve the upper limits on the velocity- weighted annihilation cross section, (av) / S 
as a function of the DM particle mass. Additionally, every time a DM particle annihilates into 
charged particles, the electromagnetic radiative correction to the main annihilation channel 
can give a more or less significant enhancement to the expected gamma-ray flux in the ob- 
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serve d environment due to internal Bremsstrahlung (IB) (lBergstromlll989l ; iBringmann et al 



20081 ). Restraining the MSSM parameters space to the stau co- annihilation region of the min- 



imal supergravity (mSUGRA) models, for instance, the wino anni hilation spectrum woul d 
receive a considerable contribution from Internal Bremsstrahlung (IBringmann et al\\200w . 
Fig. |2] shows the 95% C.L. upper limits on {av)/S as a function of the DM mass m for 
current IACTs. The projected upper limit is shown for the NFW profile, 50 h observation 
time and AQ = 2 x 10~ 5 sr. The effect of the IB is only significant below ~ 1 TeV. Some 
specific wino masses can be excluded due to the resonant enhancement in the Sommerfeld 
effect. Outside resonances, the projected upper limits are improved by more than one order 
of magnitude for DM masses above 1 TeV. The sensitivity at 95% C.L. for CTA on (o~v)/S 
as a function of the DM mass m is presented in Fig. HJ The limits are calculated for the 
ISO DM halo profile, with 200 h observation time and AQ = 10 -3 sr. The values of (o~v) 
corresponding to cosmological thermally-produced DM, (o~v) ~ 3 x 10~ 26 cm 3 s _1 , can be 
tested for specific wino masses in the resonance regions of the Sommerfeld effect. Outside 
the resonances the sensitivity on (o~v) /S is improved by more than one order of magnitude 
for TeV DM masses, reaching the level of 10~ 26 cm 3 s _1 . 



5. Astrophysical background emission 



Dwarf galaxies are generally believed to contain very little background emission from 
conventional astrophysical sources at VHE energies, and are therefore easy targets for DM 
searches. This assumption is based on their low gas content and stellar formation rate. 
However, some gamma-ray emitting sources may still exist within them: in particular from 
pulsars, and black hole accretion and/or jet emission processes. The Sagittarius and Carina 
dwarf galaxies both host globular clusters (the M54 globular cluster is located at the center of 
SgrDw), and globular clusters are known to host millisecond pulsars (MSPs). The collective 
emission of high energy gamm a-rays by MSPs in globular clusters has been detected by 
Fermi- L AT (lAbdo et aZ.ll2010bf ) . and emission in the VHE energy range has been predicted 
by several models for these objects, but has not yet been observed. The possible emission of 
very high energy radiation by millisecond pulsars from the M54 globula r cluster is examined 
in section 15.1) Add itionally, it has been suggested by some authors (see lLanzoni et al.ll2007l ; 



Noyola et al. 



20081 and references thereby) that globular clusters may host black holes with 



masses of around 10 2 to 1 4 solar masses (called intermediate-mass black holes, or IMBHs). 
Indeed, llbata et all (120091 ) suggest SgrDw may also be a possible host for a 10 4 Mq IMBH 
. Their claim is based on the study of the density profile around the central point and the 
observed rise in the velocity dispersion of stars. The high energy emission from the IMBH 
candidate in the center of M54 is discussed in section 15.21 
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5.1. Millisecond pulsars in M54 



The M54 globular cluster at the center of SgrDw is likely to harbor a large population 
of pulsars, especially MSPs. The number of MSP s in globular clusters has been shown by 
the Fermi-LAT collaboration fjAbdo et al\ 1201 Obi ) to be correlated with the collision rate V 
defined by 

r = p^rl . (9) 



In this equation, p is the central luminosity and r c is the core radi us. Taking a cen t ral su rface 
brightness of fiy ~ (14.12 — 14.9) magarcsec -2 from Table 4 of iBellazzini et all (120081 ) and 
a core radius r c = 0.9 pc , the collision rate is found to be 

r M54 ~ (o.8 - 2.6) x r M62 , (io) 

where Tm62 = 6.5 x 10 6 Lq 2 pc -2 ' 5 is the reference collision rate of the M62 glob ular cluster. 
The pr edicted number Amsp of MSPs in M54 is estimated from the collision rate (lAbdo et al 
2010bh by the relation 

Tm54 



AW = 18 + 50 x 



M62 



11) 



The collision rate from Eq. ( fit)]) gives the estimated number of MSPs in M54: Amsp = 
60 — 140. Note however that no MSP has been discovered to date in M54. 

The collective very- high-energy gamma-ray emission of millisecond pulsars from glob- 



(Bcdnarek and Sitarek 


2007D, Venter, deJag 


Cheng et al. (CCDHK) ( 


Chene et al. 


2010) 



section [2TTI one expects to observe respectively 1285 and 181 gamma-rays per hour towards 
the 47 Tucanae globular cluster, with the BS and CCDHK models. In the latter model, the 
relic gamma-rays are assumed to be the target population. The prediction of the VJC model 
is somewhat smaller, only 71 gamma-rays per hour are predicted assuming an interstellar 
magnetic field of 10 pG. The VJC model also predicts a synchrotron radiation emission. The 
emission in the keV range is predicted to be at the level of 10~ 16 TeV cm -2 s _1 for a magnetic 
field of 10 pG. This is easily accom modated by the measured d iffuse X-ray emission in M54 
which is ~ 2 x 10~ 14 TeV cm" 2 s" 1 faogdan and Gilfanovil2~o7o[ ). 



As suggested by Venter and de Jagerl (120081 ) . a rough estimate of the collective VHE 
emission of M54 can be obtained from their predicted emission of 47 Tucanae by scaling by 
the factor: 

Amsp A / ^47Tuc\ / < ^M54 > 

100 J V d M5 4 / \< W47Tuc > 



X 



(12) 
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In this equation, gL^tuc and g?m54 are the distances to 47 Tucanae and M54, and < UM54 > 
and < M47TUC > the average luminosity per cubic parsec of the globular cluster . Tak ing the 



distances, luminosity and half-mass radii of M54 and 47 Tucanae from lHarrisI (119961 ) (2010 
edition), one finds a correction factor x ~ 1.6 x 10~ 2 , assuming that M54 contains 100 MSPs. 
The expected number of gamma-rays per hour are thus 19.9 and 5.6 in the BS and CCDHK 
models. For the latter model, x was multiplied by an additional factor of 2 to take into 
account the different number of MSPs in 47 Tucanae and M54. For the VJC model, the 
number of expected gamma-rays per hour is about 1.1. 

Whether this signal is observable or not depends crucially on its spatial extension. 
The half-mass radius of M54 has an angular size of less than 1' so that the signal would 
appear almost point-like in the BS and VJC models. The CCDHK predicts an extended 
signal. The electrons and positrons responsible for the inverse Compton scattering on the 
CMB radiation have a typical diffusion length of 100 pc, which corresponds to ~ 12' at 
the distance of M54. The signal integration regions are taken as 3' for the BS and VJC 
models and 12' for the CCDHK model. With an hadron rejection factor of 10% as in section 
12.21 the number of background per hour is ~10 inside a 3' radius centered on M54. The 
significance of the collective MSP signal depends thus on the observation time T ^ s (in hours) 
as respectively 4.5 y/T ohs , 0.31 \/T obs and 0.25 y/T ohs in the BS, CCDHK and VJC models. 
The BS model would give a signal at the 4.5 a level after just a one hour observation. The 
other models would give a much smaller signal, with a typical significance of 4<r after 200 
hours of observation. 

In summary, the millisecond pulsars of M54 could give a significant VHE gamma-ray 
signal in CTA with observation times of typically 200 hours. For a cosmological thermally 
produced DM particle, (av) = 3xl0~ 26 cm 3 s~ 1 , the corresponding signal would have a 
significance of 0.1 a, after 200 hours of observation and without any boost factor. The 
collective MSP signal would be a few orders of magnitude stronger than the DM annihilation 
signal. 



5.2. Intermediate-Mass Black hole 

Significant radio and X-ray emissions are expected if the hypothesis of a central IMBH is 
valid. Unfortunately, only upper limits to the radio emission could be established from VLA 
and MOST observations. Nevertheless, these limits can be used to constrain the candidate 



black hole mass. As regards X-ray emission, Ramsay and Wu ( iRamsay and Wul l2006al Jbl) 
have analyzed the data taken with Chandra satellite. They found 7 bright sources within 
the half-mass radius of M 54. Their source number 2 lies within 1" of the density center of 
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M54. Taking into account the Chandra astrometric accuracy of 0".6 and the systematics of 
0."3 in the absolute position of t he Chandra ACS camera, this source could be associated 
with the stellar cusp identified in llbata et al\ (120091) . The source number 2 has an irregular 
shape and a luminosity of Lx = 0.72 x 10 33 ergs' - 1 faamsav and wJl2006aH ). 



The llbata et all (120091 ) estimate o f the black hole mass is consistent with the (Massive 
Black Hole - host galaxy) correlation of lFerrarese et all (120061 ) only if the host system is M54 
(M IMB hA^M54 ~ 5%). For a similar mass ratio with SgrDw as the host system, a 1000 times 
more massive black hole would be necessary, suggesting this may be a system composed of a 
dwarf galaxy hosting a prominent stellar nucleus, itself hosting a central IMBH. We estimate 
the largest contribution of the IMBH to a possible VHE gamma-ray signal, and assume that 
the IMBH is active and has a jet inclined towards the line of sight with an angle 9. The 
contribution of the black hole to th e VHE gamma-ray emission is estimated using the model 
developed by lReynoso et all (120111 ). on the emission of relativistic jets associated with active 
ga lactic nuclei. The par ameters of the model for the central black hole and jet are described 



m iReynoso et al\ ( 120111 ). The calculation also uses the constraints from the upper limits in 
the radio band and the measured X-ray emission from the source number 2. The modeled 
gamma-ray emission is shown on Fig. The parameters used in the model are given in 
Table [2j The emission depends only weakly on the black hole mass, but strongly on the 
assumed Lorentz factor Tb and inclination 9. 

The peak on the X-ray band comes from the synchrotron of electrons while a strong con- 
tribution from the Synchrotron Self-Compton (SSC) scattering can be seen at GeV energies. 
At higher energies, in particular in the CTA energy range, the emission from pp interactions 
is dominant. However, for reasonable parameters, it is in the 10~ 18 — 10~ 17 ergcm~ 2 s -1 flux 
range — too faint to be detected by CTA. 



6. Summary 



Older publications (e.g. lEvans et q/.ll2004t lAharonian et all 120081 ) on DM searches to- 



wards SgrDw used dark matter mass profiles which lead to somewhat optimistic constraints 
on particle dark matter self-annihilation cross sections. These models were used because no 
accurate modelling of SgrDw existed at that time. Several realistic models are now published 
that loosen the existing constraints by more than one order of magnitude. The future CTA 
array will be sensitive to (av) values around a few 10~ 25 cm 3 s _1 . Some models could be 
excluded after 200 hours of observation, if boosts factors are taken into account. 



However, the very high energy emission of several astrophysical objects could give an 
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observable signal for long-enough observation times. The collective very high energy emission 
of the MSPs of the M54 globular cluster, which is predicted by several models, could be much 
stronger than a DM signal. It could be observed in just a few tens of hours with CTA. The 
candidate IMBH located at the center is not expected to give an observable signal. Under 
favorable circumstances (active black hole and jet aligned towards the line of sight), it might 
nevertheless be detectable in observations of SgrDw. 
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Table 1:: Values of the LOS'-integrated squared density averaged over the solid angle (J) 
expressed in units of 10 23 GeV 2 cm -5 , for different solid angles AQ. The values of J are 
calculated for the NFW and ISO DM halo profiles. The parameters of these profiles are 
given in the first column. 



DM halo profile 


AVI = 10~ 3 sr 


Aft = 2 x 10~ 5 sr 


Afi = 2x 10~ 6 sr 


NFW 


0.065 


0.88 


3.0 


r s = 1.3 kpc 








p s = i.i xi(r 2 M pc- 3 








ISO 
r c = 0.34 kpc 
m h = 9.5 x 1O 8 M 


0.49 


1.0 


1.0 



Table 2:: Model parameters for the IMBH candidate in M54. 



Parameter 


Value 


Mbh: black hole mass 


5 x 1O 3 M 


R g : gravitational radius 


7.38 x 10 8 cm 


Lj km ^: jet kinetic power at z 


6.28 x 10 39 erg s" 1 


qy ratio 2Lj kin) /-^Edd 


0.05 


T^(zq): bulk Lorentz factor of the jet at z 


4 


9: viewing angle 


45° 


£j: jet's half-opening angle 


5° 


q rc \. jet's content of relativistic particles 


0.05 


a: hadron-to-lepton power ratio 


1 


z : jet's launching point 


50 R g 


s: spectral index injection 


2.1 


7]: acceleration efficiency 


1. x 10~ 2 


N H : column dust density 


10 21 cm" 2 
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Fig. 1. — : 95% C.L. upper limits on the velocity- weighted annihilation cross section (av) 
versus the DM mass m for a NFW (solid line) and Isothermal (ISO) (dashed line) DM halo 
profiles respectively for 11 h observation time and AQ = 2 x 10~ 5 sr. The projected upper 
limits are displayed for 50 h observation time. The sensitivities at 95% C.L. for 50 h are also 
shown for NFW (long-dashed dotted line) and ISO (dashed dotted line) DM halo profiles. 




m (TeV) 

Fig. 2. — : Projected upper limits at 95% C.L. on the (av ) /S versus the DM mass m enhanced 
by the IB (dashed line) and SE (solid line) for the NFW profile. The projected upper limits 
are shown for 50 h observation times and AQ = 2 x 10 -5 sr. 
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Fig. 3. — : Sensitivity at 95% C.L. for CTA on the velocity-weighted annihilation cross 
section (av) versus the DM mass m for a NFW (solid line) and Isothermal (ISO) (dashed 
line) DM halo profiles, respectively. The sensitivity is shown for 50 and 200 h observation 
times. The solid angle of observation is taken as Af2 = 2 x 10~ 6 sr for the NFW DM halo 
profile and AQ = 10~ 3 sr for the ISO DM halo profile. 
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m (TeV) 

Fig. 4. — : Sensitivity at 95% C.L. for CTA on the (av)/S versus the DM mass m enhanced 
by the SE for the ISO profile. The sensitivity is shown for 200 h observation times and 
Afl = 10~ 3 sr. 
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Fig. 5. — : Modeled emission of the candidate IMBH in M5 4. The inverte d empty triangles 
show the X ray emission from source 2 of Ramsay and Wu (l2006al ; l2006bl ) and the inverted 
filled triangle show radio upper limits. The various contributions to the emission are shown. 
Only the pp emission contributes in the CTA energy range. The values of the parameters of 
the black hole model are displayed in Table [2j 
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